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lit  f Abstract 

I This  report  considers  a feedback  control  law  for  linear  time-varyinp  and 

time  invariant  discrete  systems  based  on  a receding  horizon  concept  applied  to 
* a minimum  energy  problem  with  fixed  terminal  constraints.  The  control  law  is 

shown  to  be  asymptotically  stable  and  to  result  in  a new  method  for  stabilizing 
I linear  time-varying  systems  as  well  as  extending  some  well  known  methods  for 

^ stabilizing  time  invariant  systems.  In  particular,  the  stabilizing  gains  of 

1 

the  feedback  control  law  are  obtained  from  the  solution  to  a discrete  Riccati 
I equation  over  an  arbitrary  finite  time  interval,  which  is  relatively  easy  to 

compute.  The  gain  matrix  reduces  to  a constant  matrix  for  linear  time  invariant 

systems.  Some  stability  results  in  [2]  and  [5]  will  turn  out  to  be  special  cases 

of  these  results.  The  results  parallel  those  of  [4]  for  linear  continuous  time 
1 systems,  although  the  technical  details  are  tedious  and  more  involved. 


I.  Introduction 


Consider  a linear  discrete  system 

x(i+l)  = X (i)  + B(i)u(i)  (1.1) 

y(i)  = C(i)x(i)  (1.2) 

where  ,B(i)  and  C(i)  are  nxn,  nxn  and  pxn  matrices,  and  is  assximed  to 
be  nonsinpular.  Consider  also  a cost  function 


•(i)0(i)y(i)  + u'(i)R(i)u(i) 


(1.3) 


where  Q(i)  = D'(i)D(i)  i 0 and  R(i)  > 0,  together  with  the  boundary 
conditions 


x(i  ) = x 
o o 

x(i^)  = 0. 


(1.4) 


The  optimal  solution  is  obtained  by  introducing  the  2n-dimensional  Hamiltonian 
system  [11] 

n/.  \ \ f%l 


x(i+l) 

p(i) 


C’(i)Q(i)C(i), 


A 

x(i) 

p(i+l) 

(1.5) 


with  the" optimal  control 


u(i)  = -R"^(i)B’(i)p(i+l)  . 


(l.B) 


The  eouivalent  representation  of  (1.5) ^s  given  by 


» • 

x(ii-l) 

♦ j + B(i)R‘^(i)B'(i)^'‘^C'(l)Q(i)C(i),  - B(1)R‘^(1 )B*(i  )*:"^ 

x(i) 

p(ltl) 

-♦j■^C’(i)Q(i)C(l)  , 

p(l) 

(1.7) 


The  shorthand  notation  ♦j  * ♦(!♦!, i)  will  be  used  for  the  state  transition 
matrix  when  applicable. 
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Let  S(i,i^)  denote  the  2n  x 2n  state  transition  matrix  of  the  system 
(1.7)  with  the  partitioning  into  four  n x n submatrices: 


S(i,i  ) = 
o 


*(i.io) 


n(i.io) 

A(i,i^). 


The  open-loop  optimal  control  is  given  by 


(1.8) 


u(i)  = -R~^(i)B'(i) jx(i+ltiQ)-A(i+l,iQ)0~^(i^,iQ)<»(ij:,iQ^x(i^)  , (1.9) 

when  the  inverse  in  (1.9)  exists.  The  optimal  closed  loop  control  is 
given  by 

u(i)  = R’^(i)B'(i)  fl"^(i^,i+l)t(.(i^,i+l)x(i+l)  (1.10) 


whenever  the  inverse  in  (1.10)  exists.  Another  representation  of 
(1.10)  will  be  given  in  Section  111. 

The  standard  regulator  problem,  which  minimizes  the  cost  function 
(1.3)  with  a free  terminal  condition,  has  the  solution  given  by 


u(i)  = -R"^(i)B’(i)[l+K(itl,i^)B(i)R’\i)B’(i)J  K(i+l,i^)*^x(i)  . d-U) 

where  K(i,j)  is  obtained  from  ^ 

K(i,j)  = ♦j'K(i+l,j)»^  - ♦|K(itl,j)B(i)[R(i)  + B'(i)K(i+l,j)B(i^ 

X B’(i)K(i+l,j)*j  + c'(i)Q(i)C(i),  K(j,j)  = 0 . 


The  following  definition  is  necessary  for  further  analyses. 

Definition.  The  pair  {♦(i+l,i)  , B(i)}  is  said  to  be  uniformly 
completely  controllable  if  for  some  positive  integer  > 1 the  following 
conditions  hold: 

(1)  o,I<  W(i,i+1  ) < 0,1  for  all  i (1.13) 

(2)  I )| I f O3  ( |i-j  I ) for  all  i,j  * , (1.14) 

* Throughout  the  paper  the  Euclidean  norm  is  assumed  for  vectors  and  the 
spectral  norm  induced  by  the  r.uclidean  norm  for  matrices. 
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where  the  controllability  matrix  W(i,j)  is  defined  by 


W(i,j)  = I •(i,k+l)B(Jc  )B*(Ic)*'(i,k+l)  , 

k=i 


(1.15) 


and  are  positive  constants,  and  ^2^’^  maps  R into  and 

is  bounded  on  bounded  intervals. 

The  uniform  complete  observability  of  the  pair  {♦(i+l,i)  , C(i)} 
is  defined  similarly  as  above  with  the  observability  matrix 


M(i,j)  = I ♦'(k,i)C'(k)C(k)4(k,i) 
k=i 


(1.16) 


and  with  a positive  integer  I > 1 . Let  t = max  (t  ,l  } . It  is 

0 " CO 

well  known  [1,2]  that  under  uniform  complete  controllability  and  obser- 


vability conditions  together  with 


0^1  Q(i)  Ogl  and  Ogl  ^ R(i)  o^I  , 


(1.17) 


where  0^,0^,  and  are  positive  constants,  the  steady  state 

control  law  (1.11)  with  i^  = «•  is  uniformly  asymptotieally  stable,  but 
practically  speaking,  it  is  very  difficult  to  compute  K(i,«»)  for  the 
stable  control.  We  will  show  that  a modification  of  the  control  (1.10) 
results  in  an  asymptotically  stable  control  and  is  optimal  in  a certain 
sense.  In  particular,  the  gain  matrix  for  the  new  control  is  obtained  by 
solving  a Riccati  equation  on  an  arbitrary  finite  tine  interval  (larger 
than  ) , which  is  relatively  easy  to  compute  in  relation  to  the  infinite 
time  interval  for  (1.11)  with  i^  * * Preliminaries  and  some  basic 

results  are  given  in  Section  II.  The  results  for  time-varying  systems  are 
discussed  in  Section  III  and  time  invariant  systems  in  Section  IV.  A dual 
problem  will  be  discussed  in  Section  V.  Throughout  this  paper  the  following 
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matrix  identity  is  used  frequently 

I - y'(Zy'  + X)"^  Z = (I  + y'x"^Z)“^  (1.18) 

whenever  the  inverses  exist. 

II.  Some  Basic  Results. 

In  this  section  some  basic  results  are  ^iven  which  are  necessary  for 
the  succeeding  sections.  The  corresponding  results  for  continuous-time 
systems  arc  well  )cnown  FS,?  and  8].  It  appears,  however,  that  a similar 
treatment  of  discrete  systems  is  not  available  in  the  control  literature. 
Thus  we  sketch  the  proofs  briefly  in  the  Appendix. 

Theorem  2.1.  The  solution  of  the  matrix  Riccati  equation  (1.12)  is  bounded 
below  by 

- 1 ♦ tld!  ^f^  ’ (2.1) 

where  N(i,j)  is  the  solution  of  the  matrix  Lyapunov  equation 


N(i,j)  = ♦^N(i+l,j)*i  + C'(i)Q(i)C(i) 
N(j,j)  = 0 

and  the  operator  G on  , ®^^)  is  given  by 


i-1 

(Gu)(i)  = Z C(i)«(i,k+l)B(k)u(k)  . 


(2.2) 


(2.3) 


Proof.  See  Appendix  A. 

As  a consequence  of  uniform  complete  controllability,  there  exist 
positive  constants  Og  , and  such  that 

ll*(Ul,i)||  1 ag 


-1 

II  ♦ (i^l,i)||  1 »g 


|B(i)  I 


*10 


(2.4) 


I 


t 


< 


Theorem  2.2  . The  unifonn  complete  controllability  of  the  system  (1.1)  is 
invariant  under  a state  feedback  control  of  the  form 

u(i)  = K(i)x(i)  + v(i)  (2.5) 

provided  l|K(i)||  < <*^^2  some  positive  constant 

Proof:  See  Appendix  B. 

III.  Linear  Time-Varying  Systems. 

For  linear  time-varying  systems,  there  exist  few  general  methods  to 
stabilize  the  linear  system  (1.1),  one  of  which  is  the  steady  state  control 
for  (1.11)  as  mentioned  before.  We  will  suggest  another  feedback  control 
law  based  on  a recedine  horizon  concent  which  stabilized  (1.1>  and  ik  optimal 
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in  a certain  sense.  First  of  all,  we  represent  the  fixed  terminal  control 
law  (1.10)  in  terms  of  a Riccati  equation. 

Theorem  3.1.  The  fixed  terminal  optimal  closed-loop  control  law  (1.10)  can 
be  represented  as 

u(i)  = -R"^(i)B'(i)P’^(i+l,i^+l)4'.x(i)  , (3.1) 

if  the  inverse  in  (3.1)  exists,  where  P(i,j)  satisfies 
P(i,j)  = P(i+l,j)«^‘^  - P(i+l,j)4«|"^C'(i)D'(i)[I 

+ D(i)C(i)<»"^P(i+l,j)*^‘^C'(i)D'(i)]"^D(i)C(i)*:^P(i+l,i)«>:"^  (3.2) 

+ B(i-l)R"^(i-l)B’(i-l)  , P(j,j)  = 0 . 


Proof:  We  can  obtain  the  relation  (3.2)  in  a few  equivalent  ways.  It  can 
be  obtained  by  letting  K(i,j)  4 P“^(i,j)  in  Equation  (1.12)  with  K(j,j)  = “ . 
Thus  we  have 

P(i,j)  = {♦jp’^(i+l,j)^^  - ♦!p‘^(i+l,j)B(i)[R(i)  + 

B'(i)P'^(i+l,j)B(i)rt’(i)P"^(i+l,j)*.  t 
c’(i)D’(i)D(i)C(i)}‘^ 

• ♦■^{P"^(iti,j)[I-B(i)(R(i)+B’(i)P"^(i+l,j)B(i))B'(i)P'^(i+l,j)]  t 

• j‘^c'(i)D’(i)D(i)C(i)^"^}‘Hj 

• ♦■^{P"^(itl.j)[H-B(i)R"^(i)B'(i)P‘^(itl,j)]"^  + 
♦j"^c'(i)D’(i)D(i)C(i)^J^}‘^*j"^ 

X ♦'^[P(ifl.j)+B(i)R‘^(i)B'(i)]{l+*!"^c'(i)D'(i)D(i)C(i) 

• ♦■^[P(i+l,j)eB(l)R‘^(l)B’(i)]r^4j‘^  . 


iJ 
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P(i,j)  4 P(i,j)  + B(i-l)R"^(i-l)BU-l) 


(3.3) 


Then  the  above  equation  can  be  expressed  as 


P(i,j)  = ♦■^P(lel,j) 


I + c'(i)D'(i)D(i)C(i)e."^P(i+l 

1 1 


,i)}- 


+ B(i-l)R"^(i-l)B'(i-l)  . (3.4) 

Combining  (3.4)  and  (1.18)  yields  (3.2).  Since  P(i^,i^)  = 0 , P(i^,i^)  = 
B(i^-l)R”^(i^-l)Bl(i^-l)  . This  is  equivalent  to  P(i^+l,i^+l)  = 0 from 
(3.2).  Combining  (1.11)  and  (3.3)  yields  (3.1).  This  completes  the  proof. 

It  is  noted  that  the  optimal  cost  of  the  system  (1.1)  with  the 
control  (3.1)  is  given  by 


x'(i  )P"^(i  ,ip)x(i  ) . (3.5) 

o o r ^ 

The  following  lemma  is  necessary  for  the  main  theorems. 

Lemma  3.1. 

(1)  P(i,j^)  < Pd.jj)  for  i < ^3.6) 

(2)  P"^(i.jj^)  > P'^(i,j2)  for  i 1 (3-7) 

(3)  Assume  R(i)  satisfies  (1.17)  and  0 < Q(i)  < Ogl  . If  the  pair 

|♦(i+l,i)  • B(i)|  is  uniformly  completely  controllable  and  C(i)  is  bounded 
such  that  ||c(i)||  < for  all  i , then  for  a fixed  N satisfying 

t < N < • there  exist  positive  constants  and  such  that 

I lP(i,l+N)  1 Oj^gl  . (3.8) 

(4)  Assume  R(i)  and  Q(i)  satisfy  (1.17).  If  the  pairs 

{♦(i+l,i),  B(i))  and  {♦(i+l,i)  , C(i)}  are  uniformly  completely  con- 
trollable and  observable  respectively,  then  for  a fixed  N satisfying 
t < N < ••  there  exist  positive  constants  and  such  that 


I 
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Oigl  < P(i,i+N)  < . (3.9) 

Proof:  See  Appendix  C. 

From  a receding  horizon  concept  a new  control  law  is  obtained  by  re- 
placing i^  by  i+N  in  (3.1);  thus, 

u(i)  = -R~^(i)B'(i)[B(i)R"^(i)B'(i)+P(i+l,i+l+N)]"^*^x(i) 

(3.10) 

= -R~^(i)B'(i)P"^(i+l,i+l+N)*.x(i)  , N > (3.11) 

X 

A 

where  P(i,j)  may  be  obtained  from  (3.2).  Some  characteristics  of 

A ^ 

P(i,j)  are  illustrated  in  Fig.  1.  The  matrix  P(i+l,i+l+N)  is  obtained 
by  summing  (3.2)  bac)cward  fj?om  i+l+N  to  i+1  on  a finite  time  interval. 

The  most  important  property  of  the  control  law  (3.11)  is  that  it  is  a 
st^dJle  control,  though  it  is  obtained  from  a Riccati  equation  on  a finite 
time  interval. 

Theorem  3.2. 

(1) .  Assume  R(i)  satisfies  (1.17)  and  0 < Q(i)  < 0^1  . If  the 

pair  {♦(i+l,i),B(i)}  is  uniformly  completely  controllable  and  C(i)  is 
bounded  such  that  ||C(i)||  < for  all  i , then  for  a fixed  N 

satisfying  1^+1  < N < • , the  system  (1.1)  - (1.2)  with  the  feedback  control 
law  (3.11)  is  uniformly  asymptotically  stable.  (Note:  Q(i)  and  C(i) 

can  be  identically  zero). 

(2) .  Assume  R(i)  and  Q(i)  satisfy  (1.17).  If  the  pairs 
{•(i+l,i),B(i)}  and  {•(i+l,i),C(i)}  are  uniformly  completely  controllable 
and  observable  respectively,  then  for  a fixed  N satisfying 

t ♦ 1 < N < « the  system  (1.1)  -(1.2)  with  the  feedback  control  law  (3.11) 
is  uniformly  asymptotically  stable. 

Proof:  Consider  the  adjoint  system  of  (1.1)- (1.2)  with  control  law 


(3.11): 
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x(i+l)  = [♦,  - B(i)R"^(i)B*(i)P"^(i+l,i+l+N)4'^]'"^x{i)  (3.1?) 

together  with  the  associated  scalar  valued  function 

V(x,i)  = x'*TH(i+l,i+l+N)#:"^i  . (3.13) 

11 

From  Lemma  3,1  (1)  and  (2.4),  V(x,i)  satisfies 

•2  ^2  ^ 2 2 

“l4®9~  I’'!  ^ V(x,i)  5 1^1  (3.14) 

under  the  conditions  of  part  (1).  A siinilar  ineouality  involving  a,^  and 
can  be  obtained  under  the  conditions  of  part  (2).  Thus  V(x,i)  is  a 
positive  definite  function  of  x under  either  set  of  conditions.  The  dif- 
ference of  (3.13)  along  the  solution  of  the  adjoint  system  (3.12)  is  given  as 
follows : 

V(x(i),i)  - V(x(i+1) ,i+l) 

= x(i)*^^P(i+l,i+l+N)4i|  ^x(i) 

^ x'(i+l)»7^,P(i+2,i+2+N)*!“]-x(i+l) 

1+1  1+1 

= x'(i+l)[I  - B(i)R"^(i)B’(i)P"^(i+l,i+l+N)]P(i+l,i+l+N)Cl  - P'^( i+1 ,i+l+N ) 
X B(  i )r‘'^(  i )B ' ( i ) ]x(  i+1 ) - X ' ( i+i  )♦  i+2  ,i+2+N )♦  • ^^x(  i+1 ) 

= -i'(i+l)[B(i)R"^(i)B’(i)  - B(i)R‘^(i)B'(i)P"^(i+l,i+l+N)B(i)R‘^(i) 

X B'(i)]x(i+1)-  x'(i+l)[B(i)R"^(i)B'(i)  - ?(i+l,i+l+N) 

+ ♦■J;j^P(i+2,i+2+N)^|"J^]x(i+l)  . (3.15) 

From  (3.2)  we  have 

P(i+l,l+2+N)  = ♦^^j^P(i+2,i+2+N)4'["j^  + B(i)R"^(i)B'(i)  - Z(i)  , (3.16) 

where  Z(i)  is  the  non-negative  definite  second  term  on  the  right  side  of 
(3.2)  with  i and  j replaced  by  i+1  and  i+2+N  respectively. 


From  (3.16)  it  follows  that 


B(i)R"^(i)B’(i)  - P(i+l,i+l+N)  ♦ P(i+2  ,i+2+N)4>|‘J 

= P(i+l,i+2+N)  - P(i+l,i+l+N)  + Z(i)  ^ 0 , (3.17) 

where  the  last  inequality  in  (3.17)  follows  from  Lemma  3.1(1).  Thus  the 
relation  (3.15)  can  be  expressed  as 

V(x(i),i),  - V(x(i+1) ,i+l) 

S-x'(i+l)[B(i)R"^(i)B’(i)-B(i)R‘^(i)B'(i)P“^(i+l,i+l+N)B(i)p'^(i)B'(i)]x(i+l) 
= -i’(i+l)B(i)R"^''^(i)S(i)R‘^''^(i)B'(i)x(i+l)  , 


where  S(i)  is  defined  by 

_l/2  I -.1  -1/2 

S(i)AI  - R (i)B  (i)P  (i-*-l,i+l+N)B(i)R  (i)  < I . (3.18) 

If  we  can  show  that  S(i)  > for  some  positive  constant  and 

for  all  i , then  the  adjoint  system  (3.12)  will  be  asymptotically  unstable 
which  is  possible  if  and  only  if  the  original  system  is  asymptotically 
stable.  The  proof  proceeds  as  follows: 

V(x(i),i)  - V(i(itl),i+1)  < - ij^gagx'(itl)B(i)B'(i)x(i-H)  . 

Thus  we  can  have 

V(i(i+l;x  ,i  ),i+l)  - V(i  ,i  ) 
o o o o 


- “l8“8  ’‘o 


k*i 


• (i  ,ktl)B(k)B'{)c)»'(i.)c+l) 
p o’  p 


^ ®10®6“19  1^0^^  *19  * ^ *****  — ^o  ^ *c*  (3.19) 

where  ♦ (i.i  ) is  the  state  transition  matrix  of  the  closed-loop  system 
P o 
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(3.1) -(3. 11).  The  last  inequality  in  0.19)  follows  from  Theorem  2.2 
in  that  the  pair  {♦.  - B(i)R  ^(i)B'(i)P  ^(i+l,i+l+N)^j  , B(i)}  is 

uniformly  completely  controllable  since  R ^(i)B'(i)P  i+1  ,i+l+N)*^  is 

bounded  from  (1.17),  (3.8),  (3.9)  and  (2.4).  Therefore,  the  solutions  of 
the  adjoint  system  (3.12)  can  be  shown  to  increase  exponentially.  In  turn, 
this  can  be  shown  to  imply  that  the  system  (1.1)-(3.11)  is  uniformly 
asymptoticeU-ly  stable.  Thus,  it  remains  to  check  the  lower  bound  of  the 
matrix  S(i)  defined  in  (3.18).  From  (3.2)  we  have 


S(i)  = I 


r ^^^(i)B'(i)[G(i)  B(i)R'^(i)B(i)]"Vi)R‘^^^ 


I ♦ R"^^^(i)B’(i)G"^(i)B(i)R'^^^~  ^ , 


(3.20) 


where 


G(i)  = •:^,?(it2,  i+l+N){I  - ♦:"^C'(i-H)D'(i+l)[I  + Dd+DCd+D^T^  P(i+2,i  + l^N) 


i+1 


i+l 


i+1 


■•i+1^ ' ( !♦!  )D  • ( i+1 ) ] ■^D(  i+1  )C(  i+1  )*"^J^P(  i+2  ,i+l+N ) >♦ 

= P^^^( i+2 ,1+1+N)H(  i)P^^^(i+2  ,i+l+N)^' 

i+l  i+1 


(3.21) 


and  where 

A '1/2  '-1  ' ' -1  ' 

H(i)  3 I - p-^'‘‘(i+2,i+l+N)*.^JC  (i+l)D  (i+l)[I+D(i+l)C(i+l)*^^jP(i+2,i+l+N) 

K •]2l^’(i+l)D'(i+l)]'S(i+l)C(i+l)4“Jj^P^^^(i+2,  i+l+N) 


i+1 
*l/2< 


= [I  + P"-'‘‘(i+2,  i+l+N)*j"J  C’(i+l)D'(i+l)D(i+l)C(i+l)*"^^P"-''(i+2,i+l+N)]‘.^ 

(3.22) 

Note  that  G(i)  is  nonsingular  for  N ^ t^l 
From  (3.20),  (3.21)  and  (3.22)  it  follows  that 


-1 


II  S(i)||  > U+o;^ajQa-^a^^(l^a5a23o2o^5)-^]-^  S « 
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This  complatas  the  proof. ' 
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After  carefxil  inspection  of  the  matrix  Riccati  equation  (3.2)  and  the 
control  law  (3.11),  it  can  be  deduced  that  the  control  law  (3.11)  is  the 
optimal  control  law  for  the  system  (1.1)-(1.2)  which  minimizes  the  moving 
cost  fimction 

i+N-l 

I y’(k)0(>')yOc)  ♦ u'(k)R(k)u(k)  (3.23) 

k=i 

with  a moving  terminal  constraint  x(i+’N)  = 0"^  . It  will  be  interesting  to 
investigate  some  relationships  between  the  modified  control  law  (3.11)  and 
the  fixed  terminal  control  law  (3.I)  from  which  the  control  law  (3.11)  is 
obtained.  We  can  show  that  the  quadratic  cost  for  the  system  (l.l)-(1.2)  with 
the  control  law  (3.11)  is  no  more  than  that  of  the  fixed  terminal  control 
(3 .1)  with  i-  = i N. 

T O 

Theorem  3.3.  The  quadratic  cost  (1.3)  for  the  system  (1.1)-(1.2)  with  the 
control  (3.11)  has  the  following  bounds: 

x'(i  )K(i  ,i  )x(i  ) i y'(i)Q(i)ji(i)  + u'(i)R(i)u(i) 

O O i O T 
1 

o 

( x'(i  )p'^(i  ,i  ♦N)x(i  ).  (3.24) 

o 00  O 

Proof:  We  have  the  following  inequality: 

x'(i)P'^(i,l+N)x(i)  - x’(i-H)P"^(i+l,i+l+N)x(itl) 

« x’(i)P'^(i,I-H+N)x(i)  - x'(i+l)P“^(itl,i+ltN)x(i‘fl)  t 
x(i)tP’^(i,itN)  - P“^(i,i+1+N)]x(i) 
i x’(i)P"^(i,itltN)x(i)  - x'(itl)P"^(itl,i+l+N)x(i+l)  (3.25) 

Although  the  control  (3.11)  is  optimal  for  the  above  problem,  our  main  thesis 
is  that  the  receding  horizon  concept  leads  to  one  of  the  easiest  stabilization 
methods  for  linear  systems. 


I •/'Cno(£)y(i)  ♦ u'(i)R(i)u(i) 

I 

® -1  -1 

< x*(I^)P  ^(i  ,i  ♦N)x(i  ) - x’(i,)F  ♦N^xU 

O OO  O 1 li  I 

Tit*  \tMmr  bound  is  obvious.  This  completes  the  proof. 

Uhmti  the  control  law  (3.11)  is  used  as  a suboptimal  control  to  the  nti>4ilv 
etet*  control  for  (1.11),  the  error  bound  is  given  in  Theorem  3.3.  The  slab  I 
llxetlon  of  the  system  (1.1)  with  a prescribed  degree  of  stability  ran  he 
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which  follows  from  the  fact  that  P ^(i,i+N)  i P ^(i,i+l+N)  as  Riven  in 
Lemma  3.1(2).  From  (1.5)  and  (1.6)  we  can  obtain 

x'(i)p(i)  - x'(i-M)p(itl)  = x’(i)C’(i)Q(i)C(i)x(i)  + u’ (i )R(i )u(i ) (3.26) 

for  both  the  free  and  fixed  terminal  conditions  (cf.  [12]).  For  the  fixed 
terminal  condition  x(j)  = 0,  the  quantities  p(i),  u(i)  and  x(i+l)  in  (3.26) 
are  given  by 

p(i)  = P”^(i,j)x(i) 

u(i)  = -R"^(i)B'(i)[B(i)R"^(i)B'(i)+P(i+l,j)]"^«(i+l,i)x(i)  (3.27) 

x(i+l)  = [♦(i+l,i)-B(i)R'^(i)B’(i)[B(i)R"^(i)B’(i)+P(i+l,j)]'^4'(i+l,i)]x(i). 

Replacing  j by  i+l+N  in  (3.27)  and  combining  (3.26)  and  (3.25)  we  obtain 

x’(i)P"^(i,i+l+N)x(i)  - x' (i+l)P"^(i+l ,i+l+N)x(i+l) 

5 x’(i)C'(i)Q(i)C(i)x(i)  + u'(i)R(i)u(i),  (3.28) 

from  which  follows  that 

I y’(i)Q(i)y(i)  + u' (i)R(i)u(i) 
i 

® _l  _i 

( x’(i  )P  ^(i  ,i  tN)x(i  ) - x'(i,)P  ^(i,  ,i, +N)x( i , ) . 

O OO  O 1 li.  -I- 

The  lower  bound  is  obvious.  This  completes  the  proof. 

When  the  control  law  (3.11)  is  used  as  a suboptimal  control  to  the  steady 
state  control  for  (1.11),  the  error  bound  is  given  in  Theorem  3.3.  The  stabi- 
lization of  the  system  (1.1)  with  a prescribed  degree  of  stability  can  be 
obtained  from  the  result  in  Theorem  3.2. 
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Theorein  3.4.  Assume  that  the  conditions  in  Theorem  3.2  hold.  Then  the 
system  (1.1)  with  the  feedback  control  law 

u(i)  = -R"^(i)B’(i)P"^(i+l,i+l+N)f(i+l,i)x(i),  o ? 1 ( 

a 


(3.29) 


is  uniformly  asymptotically  stable,  where  P(i+l,i+l+N)  is  obtained  from  (3.2) 

a 

with  ♦(i+l,i)  replaced  by  a#(i+l,i).  Furthermore,  there  exists  an  L > 0 and 

a Y satisfying  0 < y < 1 such  that  the  transition  matrix  ♦p  (i,i^)  of  the 

a 

closed-loop  system  (l.l)-(3.29)  satisfies 


♦p  I I ^ L 


I (i-i. ) 


. i^  ? i < - 


(3.30) 


Proof:  Consider  the  system 


x(i+l)  = ot(i+l,i)x(i)  + B(i)u(i) 


(3.31) 


y(i)  = C(i)i(i)  . 


(3.32) 


It  is  easily  seen  that  the  pairs  {a'>(i+l,i),  B(i))  and  {o^(i+l,i),  C(i)} 

are  uniformly  completely  controllable  and  observable,  respectively,  if 

{♦(i+l,i),  B(i)}  and  {♦(i+l,i),  C(i)}  are  uniformly  conpletely  controllable 

and  observable  respectively.  Thus  P (i+l,i+l+N)  satisfies  the  properties  in 

a 

Lemma  3.1.  Let  (i,i  ) be  the  state  transition  matrix  of  the  system  (3.31) 
Pa  o 

with  the  control  law 


G(i)  = -R"^(i)B'(i)P'^(i+l,itl+N)a«(i+l,i)x(i)  . 

a 

Then  it  is  easy  to  see  that  ♦_  (ifi  ) - a ^o'^_  (i,i  ) . From  Theorem  3.2 

ra  O ra  o 

it  follows  that  [(ip  (i,i  )1|  J for  some  L > 0 and  0 < y < 1. 

r ^ (i-1  ) 

From  this  it  follows  that  (i,i  )||  ( lW  . This  completes  the 

o ' ^ 

proof. 


In  the  next  section  we  can  obtain  the  corresponding  results  for  linear 
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time  invariant  systems  with  simplified  forms.  Especially,  the  feedback  gain 
is  shown  to  be  a constant  matrix. 


IV.  Linear  Time  Invariant  Systems 

Consider  a linear  time  invariant  system 

x(i+l)  = ♦x(i)  + Bu(i) 
y(i)  = Cx(i) 


(4.1) 

(4.2) 


where  {♦,  B,  C}  are  constant  matrices  and  ♦ is  nonsingular,  together  with 


a cost  function 


^ y (i)Qy(i) + u (i)Ru(i) 


(4.3) 


where  0 = D D > 0 and  R > 0 are  constant  weighting  matrices.  If  the  pair 
{♦,B}  is  completely  controllable,  then  the  minimization  of  (4.3)  subject  to 
the  end  point  constraint  x(i^)  = 0 leads  to  the  optimal  feedback  control  law 


(cf.(3.1)) 


u(i)  = -R  ^B*P  ■^(i^-i)4'x(i) 


(4.4) 


if  the  inverse  in  (4.4)  exists,  where  P(k)  is  obtained  from 

P(k+1)  = ♦■^P(k)4''"^  - ♦■^P(k)*'"^C'D'(I+DC«'‘^P(k)*’"^C'D' )"^DC<»"^P(k)<f'' 
+ BR‘^B' ,P(0)  = 0. 

The  result  analogous  to  Theorem  3.2  for  this  case  is  contained  in  the 
following. 

Theorem  4.1.  Suppose  ♦ is  nonsingular.  If  the  pair  {♦.B)  is  completely 
controllable,  then  for  a fixed  integer  N satisfying  n ^ N < * the  system 
(4.1)  is  asymptotically  stable  with  the  fixed  gain  feedback  control  law 


u(i)  = -R"^B'P"^(N)»x(i) 


(4.6) 


- Ifi  - 


where  P(N)  can  be  obtained  from  (4.5)  corresponding,  to  any  chosen  pair 

{Q,R}  with  0 i 0 and  R > 0.  With  the  additional  condition  that  the 

pair  {41,0)  is  detectable  and  0 > 0,  the  result  holds  with  n f N { ■>. 

Proof:  Although  the  proof  is  merely  a specialization  of  Theorem  3.2  to  the 

time  invariant  case,  it  is  noted  that  a direct  proof  of  asymptotic  stability 

- * -1  “ -1 - 

can  be  piven  in  this  case  usinp  the  Lyapunov  function  V(x)  = x'*  P(H)*'  x 

for  the  system  x(i+l)  = [♦  - BR  ^B'P  ^{N)*]'x(i). 

The  control  law  (4.6)  is  a peneralization  of  a stabilizinp  feedback 
control  law  piven  in  [5]  involvinp  the  inverse  of  the  controllability  Gramian 
in  that  the  result  in  f5]  is  obtained  by  choosinp  Q = 0 (or  C = 0)  in  (4.5). 
That  is,  with  Q = 0,  P(N)  is  piven  by 


and 


P(N) 


I <t"^BR"  B’*’"^ 
i=0 


P‘^(N)  = ♦ 


N-1 


(4.7) 


The  advantage  of  the  control  law  (4.6)  to  the  one  employinp  (4.7)  is  that  the 
former  can  weight  the  state,  or  the  output,  by  choosinp  a proper  Q.  In  [131,  it 
is  shown  that  the  matrix  ♦ in  (4.7)  can  be  singular  for  a controllable  sinple  in- 
put system.  The  control  (4.6)  can  be  regarded  as  a general iza* ion  of  the  stable 
steady  state  control  of  (1.11)  as  shown  in  Fig.  2.  From  the  special  structure  of 
the  time  invariant  system,  the  condition  of  Theorem  4.1  can  be  weakened  as  follows 
Proposition  4.1.  Assume  t is  nonsingular.  If  the  pair  {♦,£)  is  stabilizable , 
tlien  the  system  (4.1)  is  asymptotically  stable  with  the  following  control  law 


u(i)  a -R‘^B'P^(N)»x(i),  N i n. 


(4.8) 


L7 


where  P (N)  is  the  generalized  inverse  of  the  matrix  P(N)  obtained  from 
(4.5)  for  any  Q 0 and  R > 0. 

Proof:  If  is  stabilizable , then  there  exists  a nonsinpular  real  matrix 

S such  that  with  the  transformation  x(i)  = Sx(i),  the  system  (4.1)-(4.8)  is 


transformed  to 


x(i+l)  = [♦  - BR“^B'P‘''(N)i]x(i)  , N > n. 


where  P(N)  is  obtained  from 

P(k+1)  = t”^P(k)i"^  - i"^P(k)*'"^S'~^C’D'(I+DCS‘^«'^P(k) 


)"^DCS”^i“^P(k)i'“^  + BR'^B' 


P(0)  = 0 


A A 

11  12 


0 


. B = 


and  DCS  = 


where  and  are  nonsinpular,  controllable,  and  ^ 


stable  matrix.  Direct  computation  yields 


♦ - B R"^B'P"^(N)* 

^11  “l'  “I'^ll^'^^’^ll 


x(i+l)  = 


(4.9) 


where  satisfies  (4.5)  with  ♦.DC,  and  B replaced  by  and  B^^ 

respectively.  Thus  the  matrix  in  (4.9)  is  a stable  matrix  from  Theorem  4.1. 

This  completes  the  proof. 

The  result  in  Theorem  3.4  can  be  restated  as  follows:  If  the  pair  {♦,B}  is 
completely  controllable  and  ♦ nonsinpular,  then  for  a fixed  N satisfyinp 
n { N < • the  system  (4.1)  is  asymptotically  stable  with  the  control  law 
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u(i)  = -R'^B'P"^(N)4>x(i), 
a 


a i 1 


(4.10) 


where  P (N)  is  obtained  from  (4,5)  with  ♦ replaced  by  a4>.  Furthermore,  all 
a 

poles  of  the  closed-loop  system  are  located  inside  the  circle  of  radius 
|s|  = ^ in  the  complex  plane. 

The  results  for  the  above  regulator  problem  have  other  implications  in  a 
dual  problem  of  filterinf^  theory. 

V.  Application  to  a Dual  Problem  in  Filtering  Theory. 

Consider  a linear  stochastic  system  with  white  noises, 

x(i+l)  = it'(i+l,i)x(i)  + B(i)w(i) 


y(i)  = C(i)x(i)  + v(i) 


(5.1) 


where  Ex(i  ) = x(i  ),  E(x(i  ) - x(i  ))  (x(i  ) - x(i  ))'  = T , Ew(i)  = EvC i ) = 0, 
o o o o O 0^0 

Ex(i  )w*(i  ) = Ew(i)v'(i)  = 0,  Ew(i)w'(i)  = Q(i),  and  Fv(i)v'(i)  = R( i ) . The 
o o 

standard  Kalman  filter  solution  is  given  by 


x(i+l)  = 


♦^{x(i) 


+ y(i  ,i)C’(i)[C(i)y(i  ,i)C'(i)+R(i)] 
^ o o 


-1 


[y(i)-C(i)x(i)]} 


i(i  ) = x(i  ) 
o o 


(5.2) 


where 


y(i  ,i)  is  obtained  from 
o 

y(i,j+i)  = «j(i,j){i-c'(i)cc(i)j;(i^,i)c'(i)+R(i)]’^c(i)y(i,j)}^: 


(5.3) 


+ B(i)Q(i)B*(i)  , y(i,i)  = i ? j. 

The  filtering  error  e(i)  ^ x(i)  - x(i)  has  the  mean  e'(i)  represented  by 


e(i+l)  = ♦j{I-  y(i^,i)C'(i)[C(i)  J(i^,i)C'(i)  + R(i)]'^C(i))e(i)  (5.4) 

and  variance  given  by  ^(i^.i).  The  dual  problem  of  the  fixed  terminal  minimum 

energy  problem,  (1.1)  to  (1.4),  is  considered  as  the  standard  filtering  problem 

with  a completely  unknown  initial  condition,  i.e., 

rx(i  ) = unknown  and  y = «>.  (5.5) 

o ^o 


Thus  the  filtering  problem  with  a completely  unknown  initial  condition  is  given  by 


I 

*1 

I 
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x{i+l)  = 4>j^{x(i)+r  ^(i^,i)C  (i)[C  {i)r  ^(i^,i)C(i)+R(i)]  ^ [y(i)-C(i)x(i)]  } (5.6) 

x(i  ) = arbitrary, 
o 

where 


rCi.j+l)  = ♦!"^f(i,j)*"^{I+C'(i)Q(i)C(i)»"^f(i,j)*j^}'^ 


! 

I 

I 

I 

I 

I 


r(i,j)  = r(i,j)  + c’(i)R‘^(i)C(i)  (5.7) 

r(i,i)  = 0 , 

if  the  inverse  in  (5.6)  exists.  The  error  mean  equation  for  the  estimator 
(5.6)  is  given  by 

i(i+l)  = ♦,{I-f’^(i  ,i)c'(i)R'^(i)C(i)}e(i)  . (5.8) 

1 o 

Uniform  complete  controllability  and  observability  of  the  system  (5.1) 
are  defined  as  usual  with  the  dual  system  of  (5.1).  We  state  the  corresponding 
results  in  this  case  without  proof. 

Corollary  5.1.  Assume  R(i)  satisfies  (1.18)  and  0 < 0(i)  < a^I  . If  the 
pair  {♦(i+l.i)  , C(i)}  is  uniformly  completely  observable  and  B(i)  is 
bounded , then  for  a fixed  N satisfying  + 1 < N < •*  the  state  estimator 

x(i+l)  = ♦(i+l,i)x(i)  + ♦(i+l,i)r"^(i-N,i)c'(i)[C  (i)  r~^(i-N ,i)cy i ) 

R(i)]"^[y(i)  - C(i)i(i)]  (5.9) 

is  uniformly  asymptotically  stable,  where  r(i,j)  is  obtained  from  (5.7). 

In  Corollary  5.1  the  value  of  N can  be  infinite  under  the  additional 

assumption  of  unifonn  complete  controllability  of  the  pair  {♦(i+l,i),B(i)} 

and  o I < Q(i)  < o I . It  is  easily  seen  that  the  estimator  (5.9)  is  the 
4 - “ 5 

^ I 

optimal  estimator  which  minimises  the  criterion  E(x(itl)  - x(i+l))  (x(i+l) 

- x(itl))  based  on  the  moving  information  (yd-N),  . . . y(i))  and  a 
completely  unknown  moving  initial  condition  C(x(i-M)  - x(i-tl))(x( i-N) 

f 

- x(i-N))  = * . It  is  also  noted  that  the  significance  of  the  estimator  (5.6) 


J 


?0 


lies  in  the  fact  that  it  is  one  of  the  easiest  ways  to  obtain  a linear  stable 
estimator. 


VI.  Concluding  Remarks 

An  advantage  of  the  control  law  (3.11)  is  that  the  stabilizing  feedback 
gains  are  obtained  by  summing  a Ficcati  equation  backward  in  time  over  .i 
finite  interval,  rather  than  an  infinite  time  interva.l.  The  control  law  (4.6) 
for  time  invariant  systems  generalizes  a well  known  method  of  feedback  stabili- 
zation due  to  Kleinman,  and  can  be  interpreted  as  providing  a means  for 
weighting  the  state  or  the  output  in  the  cost  function  by  choosing  Q#0  . In 

the  case  of  time  invariant  systems,  the  modified  control  law  (4.6)  can  also 
be  interpreted  as  a practical  way  to  avoid  the  singularity  near  the  terminal 
time  of  the  optimal  control  (4.4)  when  the  argument  (i^-i)  in  P ^(ij:-i) 
is  frozen  at  some  time  N = i^-i  > n . The  important  consideration  is  that 
such  expediency  still  renders  the  resulting  feedback  control  law  asymptotically 
stable.  Similar  consideration  applies  to  the  comparison  between  the  control 
laws  (3.1)  and  (3.11)  which  pertain  to  the  time  varying  system  (1-1)-(1*2). 
Corresponding  results  for  a dual  problem  in  filtering  theory  have  also  been 
discussed  in  this  paper. 


?1 


References 

[1]  R.  Z.  Kalman,  "Contributions  to  the  theory  of  --i-'toT  control,"  Bel.  "'k-. 

Mat.  Hex.,  Vol.  5/pp.  102-119,  1960. 

[2]  J.  J.  Deyst , Jr.  and  C.  F.  Price,  "Conditions  for  asymptotic  stability  of 
the  discrete  minimum-variance  linear  estimator,"  IEEE  Trans.  Automat. 

Contr. , Vol.  AC-13,  pp.  702-705,  1968. 

[3]  J.  J.  Deyst,  Jr.,  "Correction  to  the  above  paper",  IEEE  Trans.  Automat. 

Contr.,  AC-10,  pp.  562-563,  1973. 

[4]  W.  H.  l^on  and  A.  E.  Pearson,  "A  modified  minimum  enerj^y  problem  and 
feedback  stabilization  of  a linear  system".  Proceedings  of  the  Con- 
ference on  Information  Sciences  and  Systems.  Johns. Hopkins  I'niv.  Balti- 
more, MD,  March  31-April  2,  1976.  Full-lenpth  paper  submitted  for  publication. 

[5]  D.  L.  Kleinman,  "StabilizinjJ  a discrete,  constant,  linear  system  with 
application  to  iterative  methods  for  solving  the  Riccati  equation," 

IEEE  Trans.  Automat.  Contr.,  Vol.  AC-19,  pp.  252-254,  June  1974. 

[6]  B.  D.  0.  Anderson  and  J.  B.  Moore,  "New  results  in  linear  system  stability," 
SIAM  J.  Control,  Vol.  7,  No.  3,  August,  1969. 

[7]  L.  M.  Silverman  and  B.  D.  0.  Anderson,  "Controllability,  observability, 
and  stability  of  linear  systems,"  SIAM  J.  Control,  Vol.  6,  No.  1,  1968. 

[8]  R.  Canales,  "A  lower  bound  on  the  performance  of  optimal  regulators," 

IEEE  Trans.  Automat.  Contr.,  Vol.  AC-15,  No.  4,  pp.  409  - 415,  1970. 

[9]  R.  E.  Kalman  and  J.  E.  Bertram,  "Control  system  analysis  and  design  via 
the  second  method  of  Lyapunov:  II,  Discrete-time  systems,"  J.  Basic 
Engr.,  Trans.  ASME,  Ser.  D.  Vol.  82,  1960. 

[10]  P.  Dorato  and  A.  H.  Levis,  "Optimal  linear  regulator:  The  discrete-time 
case,"  IEEE  Trans.  Automat.  Contr.,  AC-16 , pp.  613-620 , DecemlM-r  , 1071 

[11]  L.  Fan  and  C.  Wang,  "The  Discrete  Maximum  Principle,"  John  Wilev  and 
Sons,  N.Y. , 1964. 


77 


References 

[12]  R.  S.  Buev,  "Two-point  boundary  value  problems  of  linear  Hamiltonian 
systems,"  SIAM  J.  Appl.  Math.,  Vol.  15,  No,  6,  November,  1967. 

[13]  W.  H.  Kwon  and  A.  E.  Pearson, "On  the  Stabilization  of  a Discrete  '"onstant 
Linear  System,"  IEEE  Trans.  Automat.  Contr.,  Vol.  AC-20,  pp,  ROO-801, 
December  1975. 


-23- 

Appendix 


A.  Proof  of  Theorem  2.1 

I 

The  proof  is  almost  the  same  as  in  [8]  for  continuous  time  systems.  Thus 


we  only  sketch  its  proof.  Vector  functions  of  time,  u(k)  and  y(k)  , 

keCi^.i^-l]  can  be  considered  as  elements  of  the  Hilbert  space  of  product 

spaces  = R™  X R™  X ..  x R™  and  x R^  x ..  x R^  (both  i^  - i^ 

times)  with  corresponding  inner  products  defined  as 

I ^f~^  , 

R " and  »y2^g  = I y^^Ck )0(k )y^(k ) . The 


k=i 


k=i 


o o 

norms  for  the  above  spaces  are  induced  by  these  inner  products.  The 


operator  G defined  in  (2.3)  maps  H.  into  H.  . Let  y.  (i)  = C(i)'Ki,i  )x 

i / h o o 

and  Jv  ■ arbitrary  control  and  y the  corresponding 

Q 

^ ^ AAA 

output.  Then  it  holds  that  y = y.  + Gu  and  J(u,x  ) = (ul^  + |Gu+y,  P . i 

h o ' ' R h p I 

The  vector  pair  (u,y)  of  the  Hilbert  space  H = H x H with  inner  product 

X £ 

aaaaaaaa  . 

<(Ui,yi)  , (u2,y2)>  = <Uj^,U2>j^  + ^yj^*y2^Q  belongs  to  the  linear  variety  j 

fs  ^ 'ft  ^ 'it  ^ 

V of  Hg  . By  the  projection  theorem,  <(u  ,y  ) , (u  - u,y  - y)>  = 0 , 

A a*  a a AVt  SV'* 

from  which  follows  that  J(u  ,x  ) = J = <(u  *u  ) , (y  ,v  ) > = < u ♦u>-, 

O r 

A * A I A|.,  . Aio 

♦ <y  ty>Q  • Also  from  the  fact  that  J = ( u P + (y  + Gu  (^  , we  can  obtain 
^ R " 0 

* I * I 7 I * I ■>  * 

J = -|Gu  I - |u  P = J + <Gu  »yf^^Q-  Prom  this  equation  follows  that 

* ^ ft  ft 

|Gu  P > (J|^  - and  f 1 - (|Gu*P/J^)(1  + 1/||g|P)  . Thus 

1 > J /J.  > 1/(1  t ||g|P)  . It  is  known  that  J = x K(i  ,i,)x  and 

“ 0“  iiii  00  x0 


= x^N(i^,i^)x^  , which  together  with  the  above  inequality  verifies  (2.1). 


- ?u  - 

B.  Proof  of  Theorem  2.2 

The  proof  is  almost  the  same  as  in  [6,7]  and  thus  a sketch  will  suffice. 
From  (2.4),  uniform  complete  controllability  implies  that  the  system  (1.1)  is 
bounded.  It  is  easily  seen  that  the  bounded  system  (1.1)  is  uniformly  com- 
pletely controllable  if,  and  only  if,  there  exists  a positive  integer  0 

such  that  W(i,i'ft^)>  for  all  i . Then  it  can  be  shown  by  contra- 

diction that  a bounded  system  is  uniformly  completely  controllable  if,  and 
only  if,  there  exists  a 1^>  0 such  that  for  every  state  ? e and  for 
any  time  i , there  exists  an  input  u defined  on  [i,i+l^-l]  such  that 
if  x(i)  = C then  x(i-»-l^)  = 0 and  |u(k)|<  y (|c|)  for  all 
kE[i,i+t^-l]  . It  is  readily  verified  that  if  the  control 
v(k)  = u(k)  - K(k)x(k)  is  the  input  to  the  feedback  closed-loop  system 
z(itl)  = [♦(!+!, i)+B(i)K(i)]2(i)  + B(i)v(i)  where  x(k)  is  the  trajectory 
of  the  open-loop  system  due  to  the  minimum  energy  control  u , then 

z<i)  = c and  z(i+l  ) = 0 (in  fact  z(k)  = x(k)  for  all  ke[i,it*,  ])  . 

c c 

Thus  we  have  |v(k)|  < |u(k)|  | I K(k)  Il|x(k ) | < Y(|r.I)  + 

k 

o,„I'Kk,i)i;  ♦ \ ♦(k,j  )B(  j )u(  j ) I < y (UI)  for  all  kc[i,i+l  -1]  where 
j = i ■ 

y(.)  and  y( . ) are  bounded  functions  on  bounded  intervals. 


The  last  inequality  of  the  above  relation  follows  from  (1.14)  and  (2.4). 
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C,  Proof  of  Lemma  3,1 


(1).  The  matrix  P(i,j)  defined  in  (3.1)  can  be  obtained  from  the 
following  free  terminal  problem: 

x(i+l)  = ♦ ^(i+l,i)x(i)  + ♦ "^(itl,i)C  (i)D  (i)u(i) 

y(i)  = B (i-l)x(i)  x(i  ) = x 

o o 

with  the  cost  function 

^ y (k)R  ^(k-l)y(k)  + u (k)u(k) 
k=i 


Since  x^P(i^,j)x^  is  the  cost  of  the  above  problem,  it  is  clear  that 
P(i,j^)  < Pd.j^)  . 

(2) .  Since  x P (i  ,j)x  is  the  minimum  cost  of  the  fixed  terminal 

o o o 

problem,  the  result  is  straightforward. 

(3) .  The  upper  bound  is  obvious  from  (1.14)  and  (2.4).  Since  P(i,j) 

comes  from  the  free  terminal  problem  given  in  (1),  we  have  from  Theorem  2.1 

P(i,itN)  > i W(i,i+N) 

■ 1 ♦ IlGlI 


where 


* -1/2  I i_'i  • I « 

(Gu)(k)  = i R (k-l)B  (k-l)»  ^(k,j)C  (j)D  (j)u(j),  kc[i,i+N-n  . 

j=l 


Thus  we  have 


-1/2 


|(Gu)(k)|  < I ||R  ' (k-l)B'(k-l)*''^(k,j)c'(j)D'(j)||  |G(j)| 

■j=i 

j=i  ® ^°[l<l<N-l  J 
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